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SUMMARY 


Methods for calculating cross sections for the breakup of galactic heavy ions by 
the Coulomb fields of the interacting nuclei are presented. With the Weizsacker- 
Williams method of virtual quanta, estimates of electromagnetic dissociation cross 
sections for a variety of reactions applicable to galactic cosmic ray shielding stud- 
ies are presented and compared with other predictions and with available experimental 
data. 


INTRODUCTION 

As the Space Station era approaches, concern is mounting over the need to pro- 
vide adequate protection for astronauts from galactic and solar cosmic rays. Al- 
though 98 percent of cosmic radiation consists of particles lighter than lithium 
(ref. 1), the relativistic nucleus component of galactic cosmic rays will be of major 
radiobiological significance for extended stays or repeated journeys into space. 

When interacting with tissue, these relativistic nuclei cause unique biological dam- 
age in the form of microlesions (ref. 2). Further, it is known that high LET (linear 
energy transfer) particles, which compose galactic cosmic rays, are highly carcino- 
genic, especially for chronic low exposures (ref. 3), and produce residual damage in 
skin many years after exposure (ref. 4). 

In previous work (refs. 5 to 17), a comprehensive nuclear interaction theory 
capable of describing absorption, total, and fragmentation cross sections at a large 
variety of energies has been developed for use as input to a radiation transport the- 
ory under concurrent development (refs. 18 to 21). This transport theory is needed 
for reliable analyses of self-shielding factors, as well as for determinations of 
personal and bulk shielding requirements. 

It has recently been found (refs. 22 to 30) that the dissociation of projectile 
nuclei by the virtual photon field of target nuclei has cross sections which are a 
sizable fraction of the nuclear projectile fragmentation cross sections. A similar 
situation also occurs for target fragmentation (ref. 29). Consequently, when compar- 
ing a theory with inclusive data, one must include, as a minimum, both the nuclear 
fragmentation process and the electromagnetic or Coulomb dissociation process. 

(These two exclusive channels may exhaust the inclusive data; although, in principle, 
one should consider other possible channels.) Thus, it is of crucial importance when 
the Coulomb dissociation cross section is a considerable fraction of the inclusive 
cross section which is true for few-nucleon removal. 

In figures 1 through 6 we have presented some simple pictures to help visualize 
the differences between dissociation due to the nuclear field and dissociation due to 
the electromagnetic field. Figure 5 shows the virtual photon field of the target nu- 
cleus interacting electromagnetically with the projectile to cause projectile excita- 
tion (and eventual breakup) . Note that this process is exactly analogous to the ex- 
citation of light nuclei induced by the electromagnetic field of an electron (fig. 7 
and ref. 31), which will be extensively studied at the 4-GeV Continuous Electron Beam 
Accelerator Facility (CEBAF) to be built in Newport News near the Langley Research 
Center. In the present investigation, the virtual photon spectrum of a target nu- 
cleus interacts with the nucleon constituents in the projectile nucleus, whereas at 



CEBAF the virtual photons of an electron will interact with the quark constituents of 
nucleons and nuclei. The energy of the virtual photons causing nuclear dissociation 
is typically on the order of 20 MeV, whereas the virtual photons at CEBAF will have 
energies up to 4 GeV (ref. 31). 

Because of the importance of nuclear electromagnetic dissociation, it is of 
great use to supplement the previously developed nuclear fragmentation theory 
(refs. 5 to 17) with calculations of the Coulomb dissociation cross section. Thus 
the present report represents an initial effort at estimating Coulomb dissociation 
cross sections. Given such a beginning effort, the methods employed here are rather 
simplistic and the resultant cross sections should be considered only as reasonable 
estimates. Specific suggestions are made as to how to improve future calculations. 

The total photodissociation cross section for removal of a particular species X 
is designated as The s y m k°l x corresponds to that defined in reference 17 

as the abladed particle in nuclear fragmentation. In general, photons, neutrons, 
deuterons, tritons, alphas, dineutrons, and so forth, will decay from a photo-excited 
nucleus; however, for the present work X is considered to be protons and neutrons 
only (i.e., one-nucleon removal). The cross section is evaluated (ref. 27) as 

°W X) = / , , V E ' X) N(E) dE (1) 

E q ( x ) 


where E Q (X) is the photonuclear threshold which actually depends on X, a v (E,X) 
is the total photonuclear reaction cross section for production of X, and N(E) is 
the virtual photon number spectrum. The calculation of N(E) and a v (E,X) is now 
considered. The symbols used in this paper are defined on pages 18 through 20. 


VIRTUAL PHOTON NUMBER SPECTRUM 

The classic Weizsacker-Williams (WW) method of virtual quanta (ref. 32) is used 
in this report. (Short discussions of this method appear in refs. 33 and 34.) 
Jackson (ref. 35) has an excellent account of this method and it is Jackson's treat- 
ment that we follow. Before proceeding, however, note that an alternative treatment 
for calculating the virtual photon spectrum of a nucleus has been presented by 
Jackie and Pilkuhn (JP) and appears in references 22, 24, and 25. The advantage of 
the JP method is that it predicts virtual photon spectra for individual multipoles, 
such as El and Ml, whereas the WW method does not. Furthermore, the JP method 
accounts for the finite extent of the charge distribution, whereas the WW method 
assumes a point charge. Olson et al. (ref. 28) provide a very clear and presentable 
discussion of the differences between the WW and the JP spectra. They note that the 
discrepancy between these two methods is not understood and must certainly be re- 
solved if further progress is to be made in this area. The minimum impact parameter 
b m in used in calculation of the virtual photon spectra is given by 


b min = R 0.1 (p) 


+ R 0>1 (T) - d 


( 2 ) 


where Rq^(P) and Rq^(t) are 1 0-percent-charge density radii of the projec- 
tile and target nuclei , * respectively (refs. 26 through 28) and d is the overlap 
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distance treated as an arbitrary parameter. Olson et al. (ref. 28) find good agree- 
ment with experimentally determined electromagnetic dissociation cross sections by 
setting d equal to 1.5 fm for the JP spectrum and to -1.5 fm for the WW spectrum. 
In fact with these values of d, one finds from table IV of reference 28 that the WW 
predictions are just as accurate, if not slightly better, than the JP prediction. 

The very similar results of these two methods is the reason for using the WW method 
in the present work. However, if one wishes to use the JP method, it is a simple 
matter to substitute the WW spectrum for the JP spectrum given on page 1531 of 
reference 28. 

The WW virtual photon number spectrum is given by 

N(E) = 1 | Z Ja ^jx K q (x) K 1 (x) - j 6 2 [kJ(x) - K*(x)]} (3 

3 

where N(E) is the number of virtual photons per unit energy E, Z t is the number 
of protons in the target nucleus, 3 is the velocity of the target in units of c, 
and a is the electromagnetic five structure constant given by 


a 



(4) 


and the parameter x in equation (3) is defined as 


x 


Eb . 


mm 


Y 3 (he) 


(5) 


where Y is the usual relativistic factor, and Kq(x) and (x) are modified 
Bessel functions of the second kind (refs. 36 and 37). The relation between the fre- 
quency spectrum dl/dE and the number spectrum is simply (ref. 35) 


N(E) 


1 n 

E dE 


(6) 


The frequency and number spectra are shown in figures 8 and 9 and are seen to be 
comparable to figure 15.8 of Jackson (ref. 35) and figure 2(a) of Olson et al. 

(ref. 28), respectively. 

As a minor technical point concerning evaluation of the Bessel functions, a gen- 
eral analytic expression for them does not exist. Jackson (ref. 35) does provide ap- 
proximate expressions for them in both the low and high frequency limits; however, in 
the present applications, these limits are not generally applicable and Jackson's ap- 
proximations fail badly. Thus, the very good polynomial approximations of Abramowitz 
and Stegun (ref. 37, pp. 378-379) are actually used here to reliably calculate the 
spectra for any frequency. 
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PHOTONUCLEAR CROSS SECTIONS 


In principle, one should really use the experimentally determined photonuclear 
reaction cross sections, as in reference 28. (Two excellent reviews of photonuclear 
reactions are given in refs. 38 and 39.) For the sake of both simplicity and gen- 
erality, however, the present work uses the parameterization of the total photo- 
absorption cross section cr abg as developed by Westfall et al. (ref. 27). The 
branching ratio g x is defined by 


VE,X) ' % » ab S (E > 


(7) 


and g x will be taken from experiment. Following Westfall et al. (ref. 27), it is 
assumed that cr abg is dominated by the electric giant dipole resonance (El GDR) 
(refs. 38 through 42) so that the present work will take a a k s to be the El GDR ab- 
sorption cross section. (This would only be approximately true (refs. 38 and 39) if 
one actually used experimental cross sections.) The absorption cross section is 
therefore given by (ref. 27) 


°abs< E) 



a 

m 


2 

E 

GDR 


) 2 /e 2 r 2 


( 8 ) 


where E GDR is the energy of the peak in the GDR cross section, T is the width of 
the El GDR, and 


a 

a = 

m irr/2 


(9) 


with the Thomas-Reiche-Kuhn cross section (ref, 43) given by 


a 

TRK 


60N t Z t 


MeV-mb 


( 10 ) 


where N fc and A^ are the target neutron and mass numbers. The GDR energy is given 
by (ref. 27) 


E GDR 


2 2 
m*c R 

o 

. 8 J 


(’ 


+ u 


he 


1 + e + 3u 
1 + £ + u 



( 11 ) 


4 



with 


u 



and 


R 


o 


r A 
o 


1/3 

t 


( 12 ) 


(13) 


where e = 0.0768, Q 1 = 17 MeV, J » 36.8 MeV, r Q = 1.18 fm, and m* is 7/10 of 
the nucleon mass. The main uncertainties in this cross section are the values of 
the branching ratios g x and the width T, which can vary from 3 to 10 MeV. The 
widths T are smallest for closed shell nuclei and largest for nonspherical nuclei 
(ref. 38). Attempts to parameterize T have not been very successful (ref. 39). 
Similarly for the branching ratios, where calculation, for instance, may involve 
knowledge of direct and statistical components as well as energy level densities of 
neighboring nuclei (refs. 40 through 42). 


Because of the uncertainties in the widths and branching ratios we have per- 
formed a detailed study by comparing theoretical cross sections with experiment as 
presented in figures 10 through 22. The aim of this comparison to experiment was to 
try to formulate an overall prescription (method) for determining F and g x which 
could be applied to systems where data do not exist. In figures 10 through 13, we 

A AC. OQ C Q 

present data and calculations for (Y,n) reactions on C, 0, Si, and Ni . The 
widths fitted were 8 MeV for 1 2 C and 10 MeV for the other three nuclei. A branching 
ratio g„ = 0.5 (suggested from equation (A7) of Westfall et al. (ref. 27)) was 

n qo i i Q7 onn 

found to be sufficient. For the (Y,n) reactions on Zr, ,ou Gd, Au, and °Pb 
(figs. 14 through 17), the widths given in figure 46 of Berman and Fultz (ref. 39) 
were sufficient as were branching ratios obtained from Weinstock and Halpern 
(ref. 44) again suggested from equation (A7) of Westfall et al. (ref. 27). For 
2 ^U(Y/n) a fit of T = 5 MeV was required (fig. 18). 


These isotopes all have a large relative abundance and we have found a general 
variation of the width with mass number also in accord with Berman and Fultz 
(ref. 39). (Note however that fig. 46 of ref. 39 is only appropriate for heavier 
nuclei which have GDR energies below a value of 18 MeV. Lighter nuclei, such as 1 6 0 
have values above this.) Thus, for naturally abundant isotopes, we feel it is safe 
to interpolate and use width values appropriate to certain mass regions as found in 
figures 10 through 18; this is done in table 1. 


The branching ratios in figures 10 through 18 are all described by what shall 
henceforth be called the branching ratio equations (BRE) defined as 
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P 


Min 


A , WH 
\ P i 


(14) 
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and 


9n = 1 - % 


where equation (14) refers to the minimum value of either Z p /A_ or the value given 
by Weinstock and Halpern in reference 44 (denoted as WH ineq. (14)). The BRE is 
suggested from equation (A7) of Westfall et al. (ref. 27) but note that their equa- 
tion is only valid in the Fe region. The BRE fits the data in figures 10 through 18 
very well. 

We warn, however, that these prescriptions for the width and branching ratios 
are not appropriate for "nonabundant" nuclei such as 1 8 0 and 54 Fe as shown in fig- 
ures 19 through 22 where both (Y,n) and (Y,p) cross sections are given. (Experimen- 
tal data for these figures are from refs. 45 and 42.) Clearly the widths are abnor- 
mal (to be suspected from ideas of the shell model) and the fitted branching ratios 
are quite different to those given by the BRE. The latter point should be obvious. 
Clearly 8 o would prefer to decay through neutron emission and ^ 4 Fe through proton 
emission. 

To summarize this comparison with the data (table 1) for "abundant" nuclei, 
values of T can be obtained from neighboring nuclei, and branching ratios can be 
determined from the BRE. For "nonabundant" nuclei, values must be obtained directly 
from experiment. 

This prescription is followed for the calculations for 20 Ne, 40 Ar, 40 Ca, 56 Fe, 
64 Cu, and 108 Ag presented in figures 23 through 28. These nuclei were chosen for 
present and future calculations of electromagnetic dissociation cross sections. 
Variations of the photoreaction cross sections with width T are shown in figures 23 
through 28. The actual values to be used in calculations are given in table 1 and 
follow the prescription of the preceding paragraph. 

Finally, in retrospect, one sees that the theory presented here for calculating 
(Y,n) and (Y,p) cross sections fits extremely well with the data (figs. 10 through 
18), given the very large mass range considered. 


ELECTROMAGNETIC DISSOCIATION CROSS SECTIONS 

As noted by Olson et al. (ref. 28), the product of the number spectrum with the 
photoreaction cross section forms a differential electromagnetic dissociation cross 
section. This cross section can be defined as 


da 


EM 


(X,E) 


dE 


= a 2 (E,X) N(E) 
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This differential cross section is finally integrated, as prescribed in equa- 
tion (1), to produce the total electromagnetic dissociation cross section. Note that 
because g x is assumed to be energy independent, we can also write 


’" ,XI = * E„ { X > 


a abs (E) N(E) dE 


g X °EM-abs (X) 


with 


a EM-abs< X > 


- I 


E o (X) 


a abs (E) 


N(E) dE 


(15) 


being the electromagnetic absorption cross section not to be confused with the photo- 
nuclear absorption cross section a a b s ^ E ^ # 

As input to the calculations, one needs the proton and neutron threshold ener- 
gies E (p) and E (n) as discussed in appendix A and listed in table 2. One also 
needs the 10-percent-charge radii discussed in appendix B and listed in table 3. The 
complete computer code listing with sample output is listed in appendix C. 

Finally one calculates the electromagnetic dissociation cross sections as listed 
in tables 4 through 7. The total (proton plus neutron) absorption cross sections for 
^Fe at 1 #88 GeV/N are given in table 4 for both d = -1 .5 and 0 fm (see eq. (2)) 
and compared with the calculations of Westfall et al. (ref. 27) who assumed d = 0. 
The reason that the present values are slightly larger than those of reference 27 is 
because they used a slightly smaller relativistic factor y to account for slowing 
down of the projectile in the target material. In table 5, comparisons are made with 
experimental values for C and 0 incident upon various targets (ref. 26). Over- 
all, one finds outstanding agreement between theory and experiment. Further, both 
values of d give comparable results. 

Unfortunately such is not the case for 1 8 0 as shown in table 6. The value of 
Up = 0.4 obtained by use of figures 19 and 20 is good but is better replaced by 

g n = 0.2 for d = -1.5 fm and by g = 0.3 for d = 0 fm. The unusual structure 
• 18 P 

in the 0(Y,n) cross section (fig. 20) may account for this discrepancy. 

1 Q7 

Target fragmentation of Au has also been studied as shown in table 7. Here 
again agreement is not as good as one would like, although the agreement is reason- 
able and better for d = 0 fm. 

In conclusion, we have been able to obtain reasonable agreement with a wide 
range of experimental results. It is suggested that a value of d = 0 fm be used in 
present and future studies. Table 8 provides a compilation of electromagnetic disso- 
ciation cross sections for use in a general fragmentation theory. Note that the 
cross section for 238 ^ on h eaV y targets is enormous. In order to improve the Coulomb 
dissociation theory, the most significant advance would be to always use the experi- 
mental photonuclear cross sections (both photoneutron and photoproton) rather than 
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calculating them as done herein. The present work has only considered neutron and 
proton removal. It would be very useful to have cross sections also for few-nucleon 
removal such as deuterons, tritons, alphas, diprotons, and dineutrons. Again experi- 
mental cross sections would be best to use. Concerning the frequency spectrum, it 
should be decided whether the WW or the JP spectrum should be used (or some other 
form) and finally the most correct value of d should be determined. 


CONCLUDING REMARKS 

Methods for calculating cross sections for the breakup of galactic heavy ions by 
the Coulomb fields of the interacting nuclei are presented. By using the Weizsacker- 
Williams method of virtual quanta, estimates of electromagnetic dissociation cross 
sections for a variety of reactions applicable to galactic cosmic ray shielding stud- 
ies are presented and compared with other predictions and with available experimental 
data. 


NASA Langley Research Center 
Hampton, VA 23665-5225 
November 15, 1985 
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APPENDIX A 


THRESHOLD ENERGIES 


For the reaction 


M p + M t * M 3 + m 4 (A1 ) 


where M p and M fc refer to the projectile and target masses, respectively, the 
projectile threshold kinetic energy for production of Mg and M^ is given by 


(M_ 


L th 


v 2 - 

2M, 


(M + 
P 


V 


(A2) 


Defining the Q-value as 


Q = (M + M t ) - (M 3 + M 4 ) 


(A3) 


equation (A2) may be written as 


L th 


-Q(M + M fc + M 3 + M 4 ) 
2M " 


(A4 ) 


For photonuclear reactions 


M p " 0 


(A5 ) 


and, therefore, for reactions like 54 Fe(y,n) and iZ S(y,&), 


32, 


T th = "2 


(A6 ) 


to a very good approximation. Note that Q is always negative for reactions because 
all reactions are endothermic, whereas decays, being exothermic, have positive values 
of Q. 
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For more bodies in the final state, such as 


M p + M t + M 3 + M 4 + M 5 + M 6 + ... + M n 


(A7) 


we simply have 


th 


(M + M + M_ + M + ... + M ) 2 - (M + \Aj 2 
3 4 5 6 N p t 

2M, 


“Q(M + M t + M 3 + M 4 + M 5 + ... + M n ) 
— 


(A8) 
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APPENDIX B 


10-PERCENT-CHARGE DENSITY RADII 

As input to the electromagnetic dissociation cross sections one requires the 
10-percent-charge density radii. De Jager et al. (ref. 46) list half-density radii 
(C) and diffuseness (z) parameters for input to density parameterizations . The 
parameter izations considered herein are the Harmonic-oscillator (HO) model , 


p(r) = p [1 + z(r/C) 2 ] exp[-(r/C) 2 ] 
o 


the 2-parameter Fermi (2pF) model. 


P(r) = 


1 + exp[ (r - C)/z] 


the 3-parameter Fermi (3pF) model. 


P(r ) = 


P q [1 + w(r 2 /C 2 )] 

1 + exp[(r - C)/z] 


and the 3-parameter Gaussian (3pG) model 


(B1 ) 


(B2 ) 


(B3 ) 


p [1 + w(r 2 /C 2 ) ] 

P (r ) = — 2 (B4) 

1 + exp(r 2 - C 2 )/z 2 


For the 2pF model one can calculate the 10-percent-charge density radius by 


R 0.1 = c + 2 « 2z (B5) 


However, such a simple analytic form is not available for the other models. Thus, 
the general method was simply to plot the various densities and determine Rq ^ 
graphically. Resultant values are listed in table 3. 



APPENDIX C 


COMPUTER CODE 


A computer program which calculates total electromagnetic dissociation cross 
sections for neutron and proton removal is given in this appendix. Required as input 
are the mass excesses of the nucleus A Z in question and also the mass excesses of 


Z and (Z - 1 ) in order to calculate proton and neutron thresholds. Further, 
the 10-percent-charge density radii, the GDR width, and the proton branching ratio 
are also required. Other inputs such as proton and mass numbers should cause no 
problem. At the end of the program is included a sample output. 
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PROGRAM LISTING 


10 

20 

30 

40 

50 

60 

70 

80 

90 

100 

110 

120 

130 

140 

150 

160 

170 

130 

190 

200 

210 

220 

230 

240 

250 

260 

270 

280 

290 

300 

310 

320 

330 

340 

350 

360 

370 

380 

390 

400 

410 

1 

420 

430 

1 

440 

1 

450 

460 

470 

430 

V) " 

490 

VI' " 

500 

510 

520 

530 

540 

550 

560 

570 

530 

590 

600 

610 


REM COULOM 

REM 

REM 

REM 

F [ XED 2 

REM 

REM 

REM NUMERICAL INTEGRATION WILL EE PERFORMED USING THE TRAPEIZG I DAL RULE 

REM 

REM 

DIM Ephot on ( 900 ) 

DIM S i grnanu ( 900 ) 

DIM Ne ( 900 ) 

REM 

REM 

REM F sc = Fine Structure Constant 

Fsc»l/137. 03604 

Hbarc=197. 32858 

Mnc <*q=938 . 95 

Mneut ron=939 . 5731 

Mprot on=938 . 2796 

Amu s =93 1.5016 

Mstar= , 7*Mncsq 

J = 36 . 8 

Q = 1 7 

Eps i 1 on=. 0768 

INPUT "ENTER GDR WIDTH (MeV)", Width 
INPUT "ENTER Z OF TARGE! T",Zt 
INPUT "ENTER A OF TARGET" , fit 
Nt =At -Zt 

INPUT "ENTER Z OF PRO JEXT I LE " , Zp 
INPUT "ENTER A OF PROJECT I LE ", Ap 
N p = A p - Z p 

INPUT "INPUT PROTON BRFiNCHING RATIO" , Frac proi. on 

INPUT "INPUT 10 percent CHARGE DENSITY RADIUS OF TARGET ( f rn > " , R 1 0t 
INPUT "INPUT 10 percent CHARGE DENSITY RADIUS OF PROJECTILE (fm)",R10p 
INPUT "INPUT Dee < overlap distance) (fm)" , Dee 
3m i n = R 1 0 1 +R 1 0p-Bee 


INPUT 

"INPUT MASS 

; EXCESS OF PROJ (MEV) 

; 

use correct 

si qri" 

, Mexcessp 

PRINT 

" < gamma, n) 

REACTION 

HAS 

NUCLEUS 

:n 

FINAL 

STATE 

WITH 

2 = ",Zp 

PRINT 

" < gamma, n ) 

REACTION 

HAS 

NUCLEUS : 

::h 

FINAL 

STATE 

WITH 

fl = " , fl p - 

PRINT 

PRINT 

" ( g am m a , p ) 

REACTION 

HAS 

NUCLEUS : 

::n 

FINAL 

STATE 

WITH 

Z = ",Zp- 

PRINT 

" ( gamma, p) 

REACTION 

HAS 

NUCLEUS : 

::n 

FINAL 

STATE 

WITH 

fl = 


PRINT 

PRINT 

PRINT 


INPUT 

, Mexc essgn 

"INPUT 

MASS EXCESS 

OF 

FINAL 

NUCLEUS 

FOR 

< g am m •=. , n ) 

REACTION 

(ME 

INPUT 

, Mexc essgp 

"INPUT 

MASS EXCESS 

OF 

FINAL 

NUCLEUS 

FOR 

< g am m s . , p ) 

REACTION 

(ME 


Mpro j = Mexc essp+Ap*Amu 
M f i n g n = M ex c e s s g n + ( A p - 1 ) * A m u 
Mf i ngp=Mexc essgp* ( Ap- 1 ) #Amu 

E t h r e =• h g n = ( ( M f i n c i n + M n e u t r o n ) A 2 - M p r o j * 2 ! 1 1 / ( 2 * M p r o j ) 

E t h r e s h g p = ( < M f i n g p + M p r o t o n > A 2 - M p r o j 2)^<2*Mproj> 

INPUT "WHAT IS KE/N OF PROJECTILE (MeV/N) ?",Tlab 
G am in a= 1 + T 1 ab / Mnc s q 
Ve 1 : =SQR (1-1 /Gamma- S 2 ) 

REM Gamma IS THE RELATIVISTIC GAMMA FACTOR OF PROJ 

REM Ve 1 IS VELOCITY OF PROJ IN UNITS OF C (RELATIVIST IC BETA FACTOR) 
Si gmam~120*Np*2p/(PI*Ap*Wi dth) 

Ro=:L . 1 :3*Ap A < 1/3) 
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620 U = 3*J*ApM-lx3>/Q 

630 Egdr=SQR < 8 . 0* J*Hbar c A 2.-’ < Mst ar*Ro A 2 ) * 1 / < l+U-< :.+Epsi 1 on+3*U>*Epsi 1 on/< 1+Epsi 

1 on+U> > ) 


REM 

REM 

REM 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 


NUMERICAL I NTEGRAT ION OR PLOT 

"neutron THRESHOLD ENERGY IS < lie V > " f Ei. hreshgn 
"prot on THRESHOLD ENERGY IS CM.?V> " , Et hreshgp 


IF Et hreshgn< Et hreshgp THEN Epho :on(l > = fit hreshgn 
II - E t h r eshgn >E t hr eshgp THEN Epho o on ( 1 > =fit hreshgp 
INPUT "ENTER ENERGY UPPER LIMIT FOR NUMERICAL INTEGRATION OR PLOT (MeVV'jE 
■'i max 

INPUT "ENTER NUMBER OF INTEGRATION OR PLOT INTERVALS" ,Npt ■= 

REM 

REM Eint is defined as. the integration or pot interval 
REM 

E i nt = < Ephot onmax-E.'phot on < 1 ) >/ CNp:- s-1 ) 

Sum = 0 

Sump=0 

Sumn=0 


R 1 = 1 TO Npts 

Ephot on=Ephot on < 1 )+< I -1 >*E i nt 
Ephoton<I)=Ephoton 

S i g m an u - S i g m am / ( 1 + ( E p h o t o n 2 - E g d r ,s 2 ') A 2 x C E p h o t o n A 2 * W i id t h A 2 ) ) 
S i grnanu < I > = S i gmanu 
E c u t o f f = H b ar c * G am m a* V e 1 / B m i n 
G = Ephot on/Ec ut of f 
CALL Bessel <G, K0, K3 > 

Ne=2*Zt *2*Fsc/< Ephot on*PI *Vel A 2>* ;G*K0*K:. 5*Ve 1 rt 2#G A 2*<Kl' 
Ne(I)=Ne 
F u n c t i o n - S i g m an u * M e 
IF 1 = 1 THEN Function =.5*Function 
IF I=Npts THEN F u n c t i o n = . 5 * F u n c t i o n 
Sum=Sum+Function 

Func t i onp = Frac prot on*Func t i on 
Func t i onn=< 1 - F'rac prot on>*Funi: t i on 
IF Ephot on< Et hreshgp THEN Fumctionp ss 0 
IF Ephot on< Et hr eshgn THEN Functionn : =0 
Sump=Sump+Func t i onp 
Sumn = Sumn + Func t i onn 


1060 NEXT I 
1070 REM 
1080 REM 
1090 REM 

1100 I nt egr al p=E i nt *Sump 

11 10 Integraln=Eint*Sumn 

1320 Integral=Integralp+Integraln 

11.30 PRINT 

1 3 40 PRINT 

13.50 PRINT "Width < Me V ) " , W i dt h 

1360 PRINT "Zt" ,Zt 

11.70 PRINT "At" ,flt 

1380 PRINT "Zp",Zp 

1190 PRINT "Ap",Ap 

1200 PRINT "KE/N < Me VxN > " , T 1 ab 

1210 PRINT "PHOTON ENERGY < NeV)" , Ephot on 

1220 PRINT 

1230 PRINT 

1240 PRINT 

1250 PRINT "Lower limit of integration < MeV Ephot on < 1 > 
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1260 

1270 

1230 

1290 

1300 

1310 

1320 

1330 

1340 

1350 

1360 

1370 

1380 

1390 

1400 

1410 

1420 

1430 

1440 

1450 

1460 

1470 

1480 

1490 

1500 

1510 

1520 

1530 

1540 

1558 

1560 

1570 

1580 

1590 

1600 

1610 

1620 

1630 

1640 

1650 

1660 

1670 

1680 

1690 

1700 

1710 

1720 

1730 

1740 

1750 

1760 

1770 

1780 

1790 

1300 

1310 

1320 

1330 

1340 

1350 

I860 

1370 

1380 

1390 

1900 

1910 


PRINT "Upper limit of integration < MeV ) " , Ephot onmax 

PRINT "Number of integration intervals is",Npts 

PRINT "Value of integration interval width (MeV)“ , Eint 

PRINT 

PRINT 

PRINT "Sigmanu < mb > " , S i gmanu 
PRINT "Sigmam < mb > " , 8 i grnam 
PRINT "Ro < f m ) " , R o 
PRINT "U",U 

PRINT " GDR Energy <MeVi",Egdr 

PRINT 

PRINT 

PRINT 

PRINT 11 PRO J VELOCITY <~Beta f ac t or > -un i t s of c" ,Vel 
PRINT "RELATIVISTIC GANNA FACTOR OF PRGJ < Me^/N > " , Gamma 
PRINT "Ecutoff ( Me V ) " , Ec ut of f 

PRINT "10 percent charge radius of target <fm) ",R10t 

PRINT "10 percent charge radius of projectile <fm)",R10p 

PRINT "Dee", Dee 

PRINT "Bmin <fm>",Bmin 

PRINT " N ( E > Cl/MeV>" , He 

PRINT 

PRINT 

PRINT "Mass excess of projectile ( MeV ) " , Mexcessp 
PRINT "Mass excess of ( pro j - neutron > < MeV > " , Mexc essgn 

PRINT "Mass excess of < pro j - proton) <MeV> ",Mexcessgp 
PRINT 

PRINT "COULOMB D I SSOC I FIT I ON CROSS SECTION (Sgrnaww) <mb)",I 
PRINT 

PRINT "Si gma< gamma, p ) < mb) ", Integra I p 

PR 1 NT " S i gma< gamma, n ) C mb > " , I nt eqr a I n 
STOP 
END 

SUB Bessel <G,K0,K1) 


Al==3. 5156229 
A2 : =3. 0899424 
A 3 ; = 1 . 2067492 
A4 : =. 2659732 
A 5 : = . 0360768 
A6 : =. 0045813 
A7 : =. 39894228 
A8 : = . 01328592 
A9==. 00225319 
A 1 0 : = . 00 1 57565 
All ; =. 00916281 
A12==. 02057706 
fil3-. 02635537 
A1 4~ . 01647633 
A15 : =. 00392377 
A16 ! = . 87890594 
fil7». 51498869 
A 1 8 : = . 15084934 
A19”. 02658733 
A20 : =. 00301532 
A21==. 0003241 1 
A22==. 39894228 
A23==. 03988024 
A24 : = . 003 6201 8 
A 25 : =. 00163801 
fi2€«. 01031555 
A27 : = . 02282967 
A2S ; = . 02895312 
A 29=. 01787654 
A 30=. 00420059 
Bl : =. 57721566 
B2==. 42278420 


n t e g r a 1 
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1920 
1930 
1940 
1950 
I960 
1970 
1980 
1990 
2000 
2010 
2020 
2030 
2040 
2050 
2060 
2070 
2080 
2090 
2100 
2110 
2120 
2130 
2140 
2150 
2160 
2170 
2180 
2190 

R13/T A 6“R:t4^T A 7+fll5/T A 8> 

2200 IF G< = 3 . 75 THEN 1 1 =G* < . 5+ft 1 6*T A 2 + fll 7*T A 4+fll 8-:fT A 6+fi 1 9*T A 3+P 20*T A 1 0 + fl2 1 *T A 1 2 

> 

2210 IF G > 3 . 75 THEN 1 1 = 1 /SQF: < G > *EXP <3 G > * < fl22-R23/T-fl24/T A 2 + fl25/7 A 3-R26/T A 4 + R27/T 
A 5-fl28/T A 6+fl29/T A 7-R30/T A 8) 

2220 S = G 2 

2230 IF G< -2 THEN K0 = -LOG < S J * I 0-B 1 + B2*S A 2 + 333*S A 4 + ]J4*S A 6 + B5*S A 84 B6*S A 1 0+B7*S A 1 2 
2240 IF G>2 THEN K0 = 1 /SQR < G ^ *EXP C -G > * < B8-B9/S+B 1 0,'S A 2-B 1 1 /S A 3 + E 1 2/SM-E 1 3/S A 5 + E 
14/S A 6> 

2250 IF G < « 2 THEN K 1 =LOG ( S ) *• 1 1 + 1 /G* < 1 +E 1 5*S A 2-B 1 6*S A 4-B 1 7*8 A 6-I 1 8*S A 8-B 19*S A 1 0- 
B 2 0 * S A 1 2 ) 

2260 IF G >2 THEN K 1 = 1 -'SQR < G 3 *EXP < -G 3> * < B2 1 t322/S-B23/S A 2 + B24/S A 3 “B25/’S A 4 + B26/S A 5 
-B27^S A 6) 

2270 SUBEND 


Note: The large array of numbers listed in the subroutine are parameters for 

determining the Bessel functions as given in reference 37. 


B3==. 23069756 
B4==. 0348859 
B5 : =. 00262698 
£6==. 0001075 
B7==. 0000074 
B8-1. 23331414 
B9«. 07832358 
j310==. 02189568 
331 1==. 01062446 
3312==. 00587872 
3313:=. 00251540 
3314==. 00053208 
3315==. 15443144 
3316:=. 67278579 
3317==. 18156897 

3318 = . 01919402 

3319 = . 001 10404 

3320 = . 00004686 
B21-1. 25331414 
3322 = . 23498619 
3323==. 03655620 
3324:=. £1504268 
3325==. 00780353 
3326:=. 00325614 
3327==. 00068245 

T=G/3.75 

IF G< =3 . 75 THEN 1 0= 1 + R 1 *T A 2+R2*T A 4+R3*T A 6+R4^T A 8+R5*T A 1 0 + F 6*T A J 2 

IF G >3 . 75 THEN I0=l/$QF:<G>*EXP<G)*<R7 + R8/T + R9/T A 2-R10/T A 34fil l/T A 4rRl 2/T A 5 + 
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SAMPLE OUTPUT 


Width <Me7> 

Zt 

At 


Zp 

Ap 

KEVN < MeV^N) 

PHOTON ENERGY <MeV> 


5.00 

02.00 
208.00 
26.00 
56.80 
1880. 00 
50.00 


Lower limit of integration <NeV) 9.67 

Upper limit of integration <NeV> 50,00 

Number of integration intervals is 108.00 

Value of integration interval width <MeV> 


.41 


Sigmanu <mb) 1.40 

Sigmam <mb) 106.41 

Ro < f m ) 4.51 

U 1.70 

GI)R Energv <MeV) 18.40 


PF!0J VELOCITY <=£eta fact or ) -uni t s of c .94 


REXATIVISTIC GRMMR 

FACTOR OF 

PRO J < MeV/N ) 

3. 80 

Ecutoff CMeV) 

42.61 



10 percent charge 

radius of 

target <fm) 

7.8 3 

10 percent charge 

radius of 

project i 1 e <f m ) 

5.28 

Dee 

0.00 



Bruin <fm) 

13.11 



NtE) Cl/MeV) 

.08 




Mass excess of projectile (NeV) 

Mass excess of <proj - neutron) <MeV> 
Mass excess of Cproj - proton) <MeV) 


-60 . 60 
-57, 4S 
-57,71 


COULOMB DISSOCIATION CROSS SEXTION <Sigmaww> <mb> 


Si gma( gamma, p) <mb) 
Si gma< gamma, n) <mb) 


258.37 

644.00 


902.37 
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SYMBOLS 


A 


BRE 


b min 

C 


c 

a 


E 

E 0 (X) 

e gdr 


GDR 




ft 


J 

K 0' K 1 

M 


m* 


N(E) 


N„ 


N + 


n 


P 


nucleon number 

nucleon number of projectile 

nucleon number of target 

branching ratio equations 

minimum impact parameter, fm 

half-density radius, fm 

speed of light, 3 x 10 8 m/sec 

overlap distance, fm 

energy, MeV 

threshold energy, MeV 

giant dipole resonance energy, MeV 

electronic charge, 1 .6 x icT 19 coul 

giant dipole resonance 

neutron branching ratio 

proton branching ratio 

branching ratio 

Planck’s constant, 6.58 x io"“ 22 MeV-sec 
intensity 

nuclear liquid drop parameter, 36.8 MeV 
modified Bessel functions of second kind 
mass, MeV/c 2 

7/10 nucleon mass, 657 MeV/c 2 

virtual photon number spectrum, MeV^ 1 

neutron number of projectile 

neutron number of target 

neutron 

projectile 
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P' 

p 

Q 

Q' 

q 

R o 

R o.i 

r 

r o 

T 

T' 

T th 

u 

V 

w 

X 

X 

Z 

z p 

z t 

z 

a 

6 

r 

Y 
e 

P(r) 


prefragment 

proton 

Q- value, MeV 

nuclear liquid drop parameter, 17 MeV 
momentum transfer 
nuclear radius, r 0 A^^ r fm 
10-percent-charge density radius, fm 
distance, fm 

radius parameter, 1,18 fm 
target 

excited target 
threshold energy, MeV 
nuclear liquid drop parameter 
speed 

nuclear density parameter 
abladed particles 
energy parameter 
proton number 

proton number of projectile 
proton number of target 
diffuseness, fm 

electromagnetic fine structure constant, 1/137 
velocity in units of c 
GDR width, MeV 
relativistic factor 

nuclear liquid drop parameter, 0.0768 

_ o 

nuclear density, fm 

— 'X 

nuclear central density, fm 
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°abs 

°EM 

°EM-abs 

a„ 

m 

0 TRK 


cross section, mb 
absorption cross section, mb 

electromagnetic dissociation cross section, mb 
electromagnetic absorption cross section, mb 
cross section parameter, mb 
Thomas -Re iche-Kuhn cross section, mb 
photonuclear cross section, mb 
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TABLE 1.- RESONANCE WIDTHS AND PARTICLE BRANCHING RATIOS 

Numbers to left of column have been confronted with experiment, numbers 
to right are our estimates used in present calculations 


Nucleus 

r, MeV 

g p 


7 Li 








Be 




12 c 

a 8.0 

b 0 .5 

b 0.5 

16 0 

a 10.0 

b 0 .5 

b 0 .5 

18 0 

a 1 2 .0 

a 0.4 

a 0 .6 

20 Ne 

c 10.0 

b 0 .5 

b 0 .5 

28 Si 

a io.o 

b 0 .5 

b 0 .5 

32 g 




40 Ar 

C 10.0 

b 0 .45 

b 0 .55 

40 Ca 

c 10.0 

b 0 .5 

b 0 .5 

48 Ti 




54 

Dq Fe 

a 3 .0 

a 0.7 

a 0 . 3 

56 Fe 

C 5 .0 

b 0 . 28 

b 0.72 

58 Ni 

a 10.0 

b 0 .5 

b 0 .5 

63 Cu 

C 5 .0 

b 0 .28 

b 0.72 

90 Zr 

d 4.0 

b 0 .05 

b 0 .95 

107 Ag 

c 5 .0 

b 0 

b 1 .0 

160 Gd 

d 4.0 

b 0 

b 1 .0 

1 81 m 




Ta 




1 ^ 7 Au 

a,d 3.5 

b o 

b 1 .0 

208 Pb 

d 3 .9 

b 0 

b 1 .0 

238y 

a 5 .0 

b o 

b 1 .0 


Fitted to data* 

D Obtained from the BRE. 
c Estimate . 

a Taken from Berman and Fultz (ref* 39) 










TABLE 2.- GIANT DIPOLE RESONANCE ENERGIES AND PARTICLE THRESHOLDS 

("Energies were calculated by equation (1); thresholds calculated"] 
L by equation (A2) J 
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TABLE 3.- THE 1 0 -PERCENT- CHARGE DENSITY RADII 


Nucleus 

10-percent radius, 
fm 

Model 

(a) 

7 Li 

3.04 

HO 

9 Be 

3.32 

HO 

1 2 C 

3.33 

HO 

16 0 

3.77 

HO and 3pF 

18 0 

3.88 

HO 

20 Ne 

4.06 

2pF 

27 Al 

4.21 

2pF 

28 Si 

4.18 

2pF 

32 S 

4.53 

3pF 

40 Ar 

4.73 

2pF and 3pF 

40 Ca 

4.80 

2pF and 3pF 

48 Ti 

5.00 

2pF 

54 Fe 

5.19 

2pF 

56 Fe 

5.28 

2pF 

58 Ni 

5.37 

2pF 

64 Cu 

5.45 

2pF 

90 Zr 

5.90 

3pG 

108 Ag and 1 07 Ag 

6.32 

2pF 

160 Gd 



181 Ta 

7.79 

2pF 

197 Au 

7.56 

2pF 

208 pb 

7.83 

2pF and 3pG 

238 0 

8.13 

2pF 


The models are defined in appendix B and are as follows 
HO, harmonic-oscillator; 2pF, 2-parameter Fermi; 3pF, 
3-parameter Fermi; and 3pG, 3-parameter Gaussian. 









TABLE 4.- CALCULATED TOTAL ELECTROMAGNETIC ABSORPTION CROSS SECTION 
FOR 1 .88 GeV/N 56 Fe INCIDENT UPON VARIOUS TARGETS 


Projectile 

Energy, GeV/N 

Target 

a EM (W) ' mb 

°EM' mb ' 

for - 




(a) 

d = -1.5 fm 

d = 0 fm 

56 Fe 

1 .88 


2 

1 .6 

1 .9 




3 

2.8 

3.3 



> 

7 

6.3 

7.3 



$ 

46 

40 

46 



63 

29 CU 

1 30 

122 

140 



107 

47 Ag 

306 





181 

„Ta 

73 

6 29 

630 

717 



208 , 
82 Pb 

834 

793 

901 



238 

92 

1008 

! 

973 

1105 


This column represents the isotope-averaged calculations of 
Westfall et al. (ref. 27). 











TABLE 5.- CALCULATED TOTAL ELECTROMAGNETIC REACTION CROSS SECTIONS FOR 
12 C AND 16 0 INCIDENT UPON VARIOUS TARGETS 



d This column represents the measurements (isotope averaged) of Heckman 
and Lindstrom (ref* 26)* 





TABLE 6.- CALCULATED TOTAL ELECTROMAGNETIC REACTION CROSS SECTIONS FOR 1 8 0 AT 

1 .7 GeV/N INCIDENT UPON VARIOUS TARGETS 


" 

Projectile 

Energy, 

GeV/N 

Target 

Final state 

<W 0) ' mb 

(a) 


o EM , mb 

, f or - 


d = 
and 

-1 .5 
g p 

fm 

of - 

d 

and 

= 0 fm 

g p of - 



0.2 

0.4 

0.3 

0.2 

18 o 

1 .7 

48 Ti 

17 0 + n 

8.7 ± 2.7 

9 

10 

12 

11 

12 

1 4 




17 n + p 

-0.5 + 1.0 

5 

4 

2 

6 

4 

3 



2°8 pb 

17 0 + n 

136 ± 2.9 

93 

108 

123 

108 

127 

144 




17 N + p 

20.2 ± 1.8 

48 

36 

24 

57 

43 

29 



238y 

17 0 + n 

140.8 ± 4.1 

113 

131 

151 

132 

154 

176 




17 N + p 

25.1 ± 1.6 

59 

44 

30 


52 

35 


a This column represents the measurements (isotope averaged) of Olson et al. 
(ref. 28). 


TABLE 7 


- TARGET FRAGMENTATION - CALCULATED TOTAL ELECTROMAGNETIC REACTION 
CROSS SECTIONS FOR VARIOUS PROJECTILES INCIDENT UPON 197 Au 


Projectile 

Energy , 
GeV/N 

Target 

Final state 

a EM (M) ' mb 
(a) 

°EM' mb ' 

for - 

d = -1.5 fm 

d = 0 fm 

12 C 

2.1 

197 Au 

198 Au + n 

66 ± 20 

33 

37 

20 Ne 

2.1 



136 ± 21 

87 

97 

40 Ar 

1 .8 



420 ± 120 

250 

278 

56 Fe 

1 .7 



680 ± 160 

488 

546 


a This column represents the data of Mercier et al. (ref. 29). 
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TABLE 8.- ELECTROMAGNETIC DISSOCIATION CROSS SECTIONS FOR A VARIETY 

OF REACTIONS WITH d = 0 fm 


Projectile 

Energy 

T, MeV 

g P 

Target 

Final state 

a EM' mb 

12 c 

86 MeV/N 

8.0 

0.5 

12 c 

"c + n 
1 1 B + P 

0.09 

0.11 


350 MeV/N 



107 Ag 

1 1 C + n 

6 






n B + p 

7 


1 .05 GeV/N 



197 Au 

1 1 C + n 
"fl + p 

31 






34 


2.1 GeV/N 



197 Au 

"c + n 
1 ’’b + p 

53 






57 

16 0 

2.1 GeV/N 

10.0 

0.5 

9 Be 

1 5 0 + n 

0.31 






15 N + p 

0.34 





12 c 

1 5 0 + n 

0.71 






15 N + p 

0.76 





2 08pb 

15 0 + n 

80 






15 N + p 

87 

40 Ar 

213 MeV/N 

10.0 

0.45 

12 C 

39 Ar + n 
39 C1 + p 

1 .2 






0.9 

56 Fe 

1 .88 GeV/N 

5.0 

0.28 

12 c 

55 Fe + n 
^Mn + p 

5.3 






2.1 





10 V 

55 Fe + n 

242 





55 Mn + p 

97 





2°8 pb 

55 Fe + n 

645 






55 Mn + p 

258 

238 0 

900 MeV/N 

5.0 

0 

12 c 

237 U + n 

237 Pa + p 

33 

0 





27 A1 

237 U + n 

142 






237 Pa + p 

0 





28 Si 

237 U + n 

165 






237 Pa + p 

0 





64 Cu 

237 U + n 

628 





4 

237 Pa + p 

0 





1 81 Ta 

237 U + n 

3208 






237 Pa + p 

0 





20Sp b 

237 U + n 

4034 






237 Pa + p 

0 















Figure 5.- Reaction diagram of projectile frag- 
mentation induced by electromagnetic 
interaction. 



Figure 6.- Reaction diagram of electromag- 
netic dissociation leading to one-nucleon 
removal. 
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o 10 20 30 40 50 


Virtual photon energy , MeV 

Figure 11.- Theoretical and experimental photoneutron reaction cross section for 
0. Width, r = 10 MeV, has been adjusted to fit data. Neutron branching 
ratio, g n = 0.5, is 1 - Z p/Ap? thus, theoretical photoneutron and photoproton 
cross sections are identical. 



10 20 30 40 50 

Virtual photon energy, MeV 


Figure 12 Theoretical and experimental photoneutron reaction cross section for 
Si. Width, r « 10, has been fitted to data. Neutron branching ratio, 

g n = 0.5, is 1 - Z /A ; thus, theoretical photoneutron and photoproton cross 
sections are identical^ 




Theoretical 

Experimental (ref. 39 ) 


Ni(Y,n) 


i' 

i \\ 


58 Ni(Y,n) = 58 Ni(y,p) 


Virtual photon energy, MeV 

Figure 13.- Theoretical and experimental photoneutron reaction cross section for 
^Ni. Width, T = 10 MeV, has been fitted to data. Neutron branching ratio, 
g n = 0.5, is 1 - Zp/Ap,- thus, theoretical photoneutron and photoproton cross 
sections are identical. 

250 - ' ■■■ Theoretical 

Experimental (ref. 39) 


Zr(y.n) 


Virtual photon energy, MeV 


Figure 14.- Theoretical and experimental photoneutron reaction cross section for 
90 Zr. Width, T = 4 MeV, has been obtained from figure 46 of reference 39. 
Neutron branching ratio, g n = 0.95, is obtained from proton branching ratio, 
g^ = 0.05, given in reference 44. 






Figure 15.- Theoretical and experimental photoneutron reaction cross section for 
Gd. Width , T = 4 MeV, has been obtained from figure 46 of reference 39. 
Neutron branching ratio, g n w 1 , is obtained from proton branching ratio, 
g p w 0, given in reference 44. 



Fi^Ygy Theoretical and experimental photoneutron reaction cross section for 

Au. Width, T = 3.5, has been obtained by fitting and from figure 46 of 
reference 39. Neutron branching ratio is g n * 1 (ref. 44)? thus, photoproton 
cross section is negligible. 
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700 



Virtual photon energy, MeV 

Figure 17.- Theoretical and experimental photoneutron reaction cross section for 
Pb. Width, T = 3.9, has been obtained from figure 46 of reference 39. 
Neutron branching ratio is g n * 1 (ref. 44); thus, photoproton cross section 
is negligible. 



Virtual photon energy, MeV 

Figure 18.- Theoretical and experimental photoneutron reaction cross section for 
U. Width, T = 5, has been fitted to data. Neutron branching ratio is 
g„ * 1 (ref. 44); thus, photoproton cross section is negligible. 



















10 20 30 40 50 


Virtual photon energy, MeV 

Figure 25.- Theoretical photoneutron and photoproton reaction cross sections 
for Ca. The branching ratios, g p - g n - 0.5 (table 1), indicate that 
photoproton and photoneutron cross sections are identical. 



0 8 16 24 32 40 


Virtual photon energy, MeV 

Figure 26.- Theoretical photoneutron and photoproton reaction cross sections for 

Fe. Width, T = 5 MeV, is taken from Westfall et al. (ref. 27) and branching 
ratios are as in table 1 . 
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Figure 27 


Figure 28, 


Virtual photon energy, lleV 

•- Theoretical photoneutron and photoproton reaction cross sections for 
4 Cu for various widths. Branching ratios are as in table 1. 



cross section is negligible; i.e., g^ * 0 (ref* 44). 
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